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ABSTRACT 

Context. There is growing observational evidence of active galactic nuclei (AGN) feedback on the ISM of radio-quiet and radio-loud 
galaxies. While AGN feedback is expected to be more common at high redshift objects, the study of local universe galaxies help to 
better characterize the different manifestations of AGN feedback. 

Aims. Molecular line observations can be used to quantify the mass and energy budget of the gas affected by AGN feedback. We 
study the emission of molecular gas in 3C 236, a Faranoff-Riley type 2 (FR II) radio source at z~0. 1, and search for the footprints of 
AGN feedback. 3C 236 shows signs of a reactivation of its AGN triggered by a recent minor merger episode. Observations have also 
previously identified an extreme H i outflow in this source. 

Methods. The IRAM Plateau de Bure interferometer (PdBI) has been used to study the distribution and kinematics of molecular gas 
in 3C 236 by imaging with high spatial resolution (0.6") the emission of the 2-1 line of 12 CO in the nucleus of the galaxy. We have 
searched for outflow signatures in the CO map. We have also derived the star-formation rate (SFR) in 3C 236 using data available 
from the literature at UV, optical and IR wavelengths, to determine the star-formation efficiency of molecular gas. 
Results. The CO emission in 3C 236 comes from a spatially resolved ~ 1 .4"(2.6 kpc)-diameter disk characterized by a regular rotating 
pattern. Within the limits imposed by the sensitivity and velocity coverage of the CO data, we do not detect any outflow signatures in 
the cold molecular gas. The disk has a cold gas mass M(H2)~2.1xl0 9 Mq. Based on CO we determine a new value for the redshift of 
the source zco=0.09927±0.0002. The similarity between the CO and Hi profiles indicates that the deep H i absorption in 3C 236 can 
be accounted for by a rotating H i structure. This restricts the evidence of H i outflow only to the most extreme velocities. In the light 
of the new redshift value, the analysis of the ionized gas kinematics reveals a fast (~1000 km s~') outflow. As for the CO emitting 
gas, outflow signatures are nevertheless absent in the warm molecular gas emission traced by infrared H2 lines. The star-formation 
efficiency in 3C 236 is consistent with the value measured in normal galaxies, which follow the canonical Kennicutt- Schmidt (KS) 
relation. This result, confirmed to hold in other young radio sources examined in this work, is in stark contrast with the factor of 10-50 
lower SFE that has been claimed to characterize evolved powerful radio galaxies. 

Conclusions. There are no signs of ongoing AGN feedback on the molecular ISM of 3C 236. The recent reactivation of the AGN in 
3C 236 (about ~10 5 yr ago) is a likely explanation for the early evolutionary status of its molecular disk. 

Key words. Galaxies: individual: 3C 236 - Galaxies: ISM - Galaxies: kinematics and dynamics - Galaxies: active - ISM: jets and 
outflows 



1. Introduction 

1.1. AGN Feedback 

Active Galactic Nuclei (AGN) release vast amounts of energy 
into the interstellar medium (ISM) of their host galaxies. This 
energy input can heat gas in the ISM, preventing its collapse and 
inhibiting star formation. It may also expel the gas of the ISM 
in the form of winds (outflows) which deplete the host of star- 
forming fuel. This transfer of energy from the AGN to the host is 



Send offprint requests to: 

Alvaro Labiano: alvaro . labiano@cab . inta-csic . es 

Based on observations carried out with the IRAM Plateau de Bure 
Interferometer. IRAM is supported by INSU/CNRS (France), MPG 
( Germany ) and IGN ( Spain ). 



known as AGN feedback. AGN feedback could be responsible of 
the correlations between black hole and host galaxy bulge mass 
(Magorrian et al. 1998; Tremaine et al. 2002; Marconi & Hunt 
2003; Haring & Rix 2004). AGN feedback can also explain the 
fast transition of early type galaxies from the blue-cloud to the 
red-sequence (Schawinski et al. 2007; Kaviraj et al. 201 1). 

Over the last decade, AGN feedback has been increasingly 
included in models of galaxy evolution (King 2003; Granato 
et al. 2004; Di Matteo et al. 2005; Croton et al. 2006; Ciotti 
& Ostriker 2007; Menci et al. 2008; King et al. 2008; Merloni 
& Heinz 2008; Narayanan et al. 2008; Silk & Nusser 2010). 
Observational searches for signatures of AGN feedback have 
also grown in number (Thomas et al. 2005; Murray et al. 2005; 
Schawinski et al. 2006; Miiller Sanchez et al. 2006; Feruglio 
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et al. 2010; Crenshaw et al. 2010; Fischer et al. 2011; Villar- 
Martfn et al. 2011; Dasyra & Combes 2011; Sturm et al. 2011; 
Maiolino et al. 2012; Aalto et al. 2012). Yet, observational evi- 
dences of AGN feedback are still fragmentary. 

AGN feedback takes place through two mechanisms: the ra- 
diative or quasar mode, and the kinetic or radio mode. In the 
quasar mode, the radiation from the AGN dominates the energy 
transfer to the ISM. In the radio mode, the momentum of the jet 
is transferred to the ISM. 

The most dramatic effects of the kinetic mode are seen in 
powerful radio galaxies. With jets up to several Mpc-size, they 
can inject enormous amounts of energy, not only in their host 
galaxy ISM, but also in the IGM of the galactic group or cluster 
where they reside (e.g., Birzan et al. 2004; Fabian et al. 2006; 
McNamara & Nulsen 2007, 2012). It is thought that a signifi- 
cant fraction of massive galaxies undergo a radio galaxy phase at 
least once in their lifetime (e.g., Best et al. 2006). Consequently, 
understanding how AGN feedback works during this phase is 
crucial to a better understanding of galaxy evolvution. 

Observations of H i and ionized gas in radio galaxies have 
found massive outflows in a significant number of radio galax- 
ies at different redshift ranges (e.g., Morganti et al. 2003a; 
Rupke et al. 2005; Holt et al. 2006; Nesvadba et al. 2006, 2008; 
Lehnert et al. 201 1). Even though molecular gas may dominate 
the mass/energy budget of the wind, unambiguous evidence of 
a molecular outflow has only been found thus far in one radio 
galaxy (4C 12.50; Dasyra & Combes 201 1). Another manifesta- 
tion of AGN feedback, the inhibition of star formation, has also 
been claimed to be at work in radio galaxies. Nesvadba et al. 
(2010) used the 7.7 /mi polycyclic aromatic hydrocarbons (PAH) 
emission to estimate the star-formation rate (SFR) in a sample of 
radio galaxies. Based on their estimated molecular gas contents, 
Nesvadba et al. (2010) found that radio galaxies in their sam- 
ple are ~ 10-50 times less efficient in forming stars compared to 
normal galaxies that follow the canonical Kennicutt-Schmidt re- 
lationship (Schmidt 1959; Kennicutt 1998). The warm-H2 emis- 
sion of these radio galaxies suggests the presence of shocks in 
the warm molecular gas traced by a set of IR H2 lines. These 
shocks, identified in the warm-H2 phase, would increase the tur- 
bulence in the molecular ISM as a whole, inhibiting the star for- 
mation in the host. However, it is still an open issue whether 
shocks are affecting the bulk of the molecular ISM or if they only 
concern the warm molecular gas phase. CO observations, well 
adapted to trace most of the molecular ISM, are key if we are 
to understand the complexity of how AGN feedback affects the 
kinematics and the star-formation properties of molecular gas in 
radio galaxies. 

1.2. 3C236 

3C 236 is the second largest radio galaxy in the Universe 
(~4.5 Mpc deprojected size, Willis et al. 1974; Barthel et al. 
1985) 1 . Its radio structure shows an old (2.6 x 10 8 yr), classi- 
cal double, FR II source (Fanaroff & Riley 1974) on large scale 
maps, and a younger (10 5 yr), ~2 kpc Compact Steep Spectrum 
source (CSS, O'Dea 1998) in the center, responsible of two 
thirds of the total radio emission of 3C 236 (Schilizzi et al. 
2001). Both structures are roughly aligned and oriented at an 
angle ~30° to the plane of the sky, with the North- West jet ap- 
proaching (Schilizzi et al. 2001). The double-double morphol- 



1 Machalski et al. (2008) recently reported the discovery of the cur- 
rently largest radio galaxy, J1420-0545, with a projected size of 4.7 
Mpc. 



ogy (CSS + large-scale source; e.g., Schoenmakers et al. 2000; 
Kaiser et al. 2000) of 3C 236 is consistent with a reignition of 
the AGN activity, probably due to a minor merger (O'Dea et al. 
2001). 

The host of 3C 236 is a massive (~10 12 ~ 13 Mq, Sandage 
1972; Strom & Willis 1980) elliptical galaxy (Barthel et al. 1985; 
Mezcua et al. 2011), with distorted optical morphology at kpc 
scales (Smith & Heckman 1989), consistent with a merger. High 
resolution imaging and optical broadband-filter absorption maps 
(de Koff et al. 2000) show an inner dust disk (0.5 kpc from the 
nucleus), and a broad, dynamically young, ~8 kpc long dust lane 
(1.5 kpc from the nucleus), misaligned by ~25° with the inner 
disk (Martel et al. 1999; de Koff et al. 2000). The total mass of 
the whole dust system is ~10 7 Mq (de Koff et al. 2000; Sodroski 
et al. 1994). 

The total stellar mass of 3C 236 has been estimated to be 
10 12 Mq (Tadhunter et al. 201 1), consisting of a young (<10 7 yr) 
and an old stellar population (>10 Gyr). The old stellar pop- 
ulation represents 72% of the total stellar mass (see also Holt 
et al. 2007; Buttiglione et al. 2009). Using high resolution opti- 
cal and UV photometry O'Dea et al. (2001) and Tremblay et al. 
(2010) found young (~10 7 yr) star forming knots on the edge 
of the dust lane, surrounding the nucleus, as well as and older 
(~10 9 yr) population in the nucleus. The star formation in the 
young knots could have been triggered by the infalling gas from 
the minor merger which reignited the AGN. However, it is not 
completely clear how the nuclear population relates to the AGN 
activity (Tremblay et al. 2010). 

The star-formation history of 3C 236 and the possibly re- 
lated reignition of its AGN make 3C 236 an ideal candidate for 
searches of AGN feedback signatures. Furthermore, 3C 236 is 
known to have a fast (~1500 km s~'), massive (~50 Mq yr -1 ) 
H 1 outflow: Westerbork Radio Telescope spectra show a shallow, 
narrow blueshfited H 1 absorption component with a broad blue 
wing (Morganti et al. 2005). Previous searches of molecular gas 
in 3C 236 were plagued by limited sensitivity and they were un- 
dertaken with single-dish telescopes. Based on CO observations 
done with the IRAM-30m telescope in a sample of radio galax- 
ies, Saripalli & Mack (2007) reported the non-detection of CO 
line emission in 3C 236; the implied upper limit on the molecular 
gas mass of 3C 236 (~2-3xl0 9 Mq) was not very compelling, 
however. 

This is the second of a series of papers where we use interfer- 
ometric observations of the molecular line emission to study the 
fueling and the feedback of activity in a sample of nearby radio 
galaxies (observations of 4C 31 .04, a young compact symmet- 
ric object (CSO) source, were presented by Garcia-Burillo et al. 
2007). To this purpose, we have carried out 1 mm/3 mm contin- 
uum and 12 CO(2 - 1) line high resolution and high sensitivity 
observations of 3C 236 with the Plateau de Bure Interferometer 
(PdBI). We study the continuum emission of the source at 1 mm 
and 3 mm wavelengths and analyze the distribution and kinemat- 
ics of molecular gas based on the emission of the 12 CO(2 - 1) 
line. We search for outflow signatures in the CO map and re- 
evaluate the evidence of outflow signatures in different ISM 
tracers. Furthermore, we derive the star-formation rate (SFR) 
of 3C 236 based on several tracers available from the litera- 
ture. Based on the gas mass derived form CO emission and on 
the SFR estimates, we compare the star-formation efficiency of 
3C 236 against a sample of powerful radio galaxies. Finally, we 
discuss a scenario consistent with the properties and the history 
of 3C 236 and its host. 
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Fig. 1. Continuum maps of 3C 236 at 1 mm (left) and 3 mm (right). Crosses (+) mark the position of the AGN and jet peak emission according to 
the UV-FIT models (Table 1). Triangles mark the position of components A and Bl from the VLBI map of Schilizzi et al. (2001). The contour 
levels are -3 mjy beam -1 , 3 mjy beam -1 , 9 mjy beam -1 to 63 mjy beam -1 in steps of 9 mjy beam -1 for the 1 mm map, and -6 mJy beam -1 , 6 
mjy beam -1 , 22 mjy beam -1 to 176 mJy beam -1 in steps of 22 mjy beam -1 for the 3 mm map. The filled ellipses represent the beams at 1 mm 
(0.60"x0.50", PA=41°) and 3 mm (1.16"x0.86", PA=41°). (As, A<J)=(0,0) corresponds to the position of the AGN at 1 mm in both panels. 



Table 1. Point source models. 



Component 


Wavelength 


RA (J2000) 


Dec (J2000) 


Flux (mJy) 


AGN 


1 mm 


10 06 01.753 


+34 54 10.425 


71.4 ±0.4 


Jet 


1 mm 


10 06 01.717 


+34 54 10.777 


10.1 ±0.4 


AGN 


3 mm 


10 06 01.748 


+34 54 10.385 


150.4 ± 0.8 


Jet 


3 mm 


10 06 01.699 


+34 54 10.758 


45.2 ± 0.8 


VLBI core 


6 cm 


10 06 01.756 


+34 54 10.460 





AGN and jet positions in the 1 mm and 3 mm maps of 3C 236, 
according to the UV-FIT models. The VLBI position of the core 
(Schilizzi et al. 2001; Taylor et al. 2001) is included for comparison. 



2. Observations 

2.1. Interferometer observations 

PdBI (Guilloteau et al. 1992) observations of 3C 236 were ob- 
tained in the 12 CO(2 - 1) 1mm line using six antennas in the 
array B's configuration in January 2009. We assumed the op- 
tical redshift of the source derived by Holt (2005), z=0. 10054 
(uo(HEL)=30 129 km s -1 ), to tune the 1 mm receivers centered 
on the 12 CO(2-l) line redshifted to 209.504GHz. The relative 
velocity scale, although re-determined in Sect. 4, initially refers 
to this redshift (v-u ). With this setting we obtained a velocity 
coverage of ~1400 kms -1 at 209.5 GHz with the narrow band 
correlator of the PdBI (1 GHz- wide) and the two polarizations 
of the receiver. Observations were conducted in a single pointing 
of size ~23". The adopted phase-tracking center of the obser- 
vations was set at (a 2 ooo, e5 2 ooo)=(10 h 06 m 01.7 s , 34°54'10"), the 
position of the nucleus given by NED. Nevertheless the position 
of the dynamical center determined in Sect. 3, which coincides 
with the position of the radiocontinuum VLBI core (Schilizzi 
et al. 2001), is ^1" offset to the NE with respect to the ar- 
ray center: (A a, A5)~(0.6", 0.4"). Visibilities were obtained 
through on-source integration times of 22.5 min framed by short 
(2 min) phase and amplitude calibrations on the nearby quasars 
0923+392 and 1040+244. The visibilities were calibrated using 
the antenna-based scheme. The absolute flux density scale was 
calibrated on MWC349 and found to be accurate to <10% at 
209.5 GHz. 



The image reconstruction was done with the standard 
IRAM/GILDAS software (Guilloteau & Lucas 2000) 2 . We used 
natural weighting and no taper to generate the CO line map 
with a size of 66" and 0.13" /pixel sampling, and obtained a 
synthesized beam of 0.60"x0.51" @PA=41° at 1 mm. The lcr 
point source sensitivity was derived from emission-free chan- 
nels resulting in 1.1 mjy beam -1 in 20MHz(~29kms -1 )-wide 
channels. We obtained a map of the continuum emission at 
209.5 GHz averaging channels free of line emission from v- 
u o =+100 to +640 kms -1 . The corresponding l<x sensitivity of 
continuum emission is ~0.3 mjy beam -1 at 209.5 GHz. 

We subsequently observed 3C 236 in the HCO + (1-0) and 
HCN(l-O) transitions at 3 mm for six hours on January 2010. 
For these observations we tuned the receivers at a frequency in- 
termediate between those of the HCO + (1-0) (89.188 GHz) and 
HCN(l-O) (88.632 GHz) transitions redshifted to 81.017GHz. 
The velocity coverage at 8 1 GHz was -3700 kms -1 . We used 
the same phase tracking center and redshift as adopted for the 
1 mm observations described above. The twofold goal of the 
3 mm observations was the detection of the HCO + (1-0) and 
HCN(l-O) lines as well as of their underlying continuum emis- 
sion. Notwithstanding, these lines were not detected down to 
the sensitivity limit (0.8 mjy beam -1 in 10 MHz, ~37 km s -1 
-wide channels). A 3 mm continuum map was built inside the 
59" field-of-view using velocity channels away from any poten- 
tial contribution of line emission in the signal side band of the 
PdBI receivers (u-u o =-1370 to -1740 kms -1 ). The corresponding 
lcr sensitivity is ~0.3mJybeam -1 at 81 GHz. We used natural 
weighting to generate the maps of the continuum emission and 
obtained a synthesized beam of 1.17"x0.86" @PA=Al°. 

2.2. Ancillary data 

We have used the following archival HST images: 
ACS/SBC/F140LP (hereafter FUV image), ACS/HRC/F555W 
(V-band), STIS/NUV-MAMA/F25SRF2 (NUV, O'Dea et al. 
2001; Allen et al. 2002; Tremblay et al. 2010); WFPC2/F702W 

2 http://www.iram.fr/IRAMFR/GILDAS 
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Fig. 2. I2 CO(2 - 1) intensity map of 3C 236, integrating all emission 
above 3tr levels, from d-d„=-685 to +165 km s . The contour levels are 
3cr, to 6.5<r in steps of 0.5cr, with cr=0.21 Jy beam -1 km s - '. (Aa, A<5)- 
offsets in arcsec are relative to the location of the AGN. The gray ellipse 
shows the 1 mm beam. 



(fl-band) and NICMOS2 (//-band) obtained by Sparks and 
collaborators (de Koff et al. 1996; McCarthy et al. 1997; Madrid 
et al. 2006; Tremblay et al. 2007; Floyd et al. 2008); Hi data 
by Morganti et al. (2005); mid-IR data from Dasyra & Combes 
(2011) and Guillard et al. (2012); Sloan Digital Sky Survey 
(SDSS) spectra from York et al. (2000), Abazajian et al. (2009), 
and Buttiglione et al. (2009). 

We use H = 71, Q M = 0.27, Q A = 0.73 (Spergel et al. 
2003) throughout the paper. Luminosity and angular distances 
are Di=452Mpc and D,4=374Mpc; the latter gives a physical 
scale of l"=1.8kpc (Wright 2006). 

3. Continuum maps 

Figure 1 shows the continuum maps derived at 209.5 GHz and 
81 GHz in the nucleus of 3C 236. The emission spreads on spa- 
tial scales significantly larger than the PdBI synthesized beam 
at both frequencies and consists of a dominant central com- 
ponent and a fainter emission component that extends to the 
NW of the source. We used the GILDAS task UV-FIT to fit 
the continuum visibilities at both frequencies with a set of two 
point sources. Table 1 lists the parameters of the best fit solu- 
tions. Most of the flux comes from a central source located at 
(Aa, A(5)~(+0.6+0.1", +0.4+0.1") at both frequencies. This cor- 
responds within the errors with the position of the AGN core 
determined in the VLBI 1 .7 and 5 GHz maps (component B2 
in Schilizzi et al. 2001; Taylor et al. 2001) and is therefore 
adopted as the best guess for the dynamical center of 3C 236 
(component labelled as 'AGN' in Fig. 1). The NW 1 mm and 
3 mm components are both fitted by a point source located at 
(Aa, A5)~(+0.0+0.1", +0.8+0.1"), i.e., at a position interme- 
diate between the jet components A (radio lobe) and Bl (base 
of the radio jet) identified in the VLBI 1.7 and 5 GHz maps of 
Schilizzi et al. (2001) and Taylor et al. (2001). This indicates that 
the NW elongations identified in Fig. 1 are the mm-counterparts 
of the approaching radio jet of the inner region of 3C 236 (com- 
ponent labelled as 'jet' in Fig. 1). Using the fitted positions for 
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Fig. 3. Spectrum of the 12 CO(2 - 1) line emission (after subtracting 
the continuum emission) towards the AGN position of 3C 236 (his- 
togram). The red solid line shows the Gaussian fit to the spectrum. 
Velocities in the x-axis are with respect to the redshift z = 0.10054 
(do(HEL)=30 129 km s~'), where the PdBI narrow band correlator was 
centered. Arrows mark the heliocentric systemic velocity of 3C 236 
according to SDSS (29 740 km s" 1 ), CO (29 761 km s -1 ), and NED 
(30 129 km s- 1 ). 



the 'AGN' and 'jet' components we derive a PA of 302° for the 
mm-jet(s); this is consistent with the orientation of the CSS radio 
jet determined in the VLBI cm-maps of the source. 

To establish the nature of the 1 mm and 3 mm emissions 
of 3C 236 we have compared the fluxes of the fitted compo- 
nents ('AGN' and 'jet') with the corresponding radio emissions 
between 10 MHz and 10 GHz of 3C 236 (available in NED). 
In the log f v - log v plot, the 1 mm and 3 mm data points 
are aligned with the radio data points, following a power law 
with spectral index ^-0.9+0.1. This spectral index is expected 
for synchrotron radiation from CSS sources (e.g., O'Dea 1998). 
Therefore, we conclude that the 1 mm and 3 mm emission in 
3C 236 can be accounted for by synchrotron radiation. 

4. CO maps 

4. 1 . Distribution of molecular gas 

We detect significant CO(2-l) emission above 3<x levels from 
i>-t; =-685 to +165 km s _1 . Figure 2 shows the CO(2-l) inten- 
sity map obtained by integrating the emission within this veloc- 
ity range with no threshold value adopted on the intensities. As 
illustrated in this figure, CO emission is detected at significant 
levels only in the inner diameter of 2" (3.6 kpc) of 3C 236. The 
distribution of the molecular gas can be described as a spatially- 
resolved elongated disk-like source. While the peak of CO emis- 
sion lies close to the AGN, the detection of lower level emission 
that extends further to the NE gives the overall impression that 
the molecular disk, which has a deconvolved major axis diame- 
ter of 1.3" (2.3 kpc), is off-centered with respect to the AGN. 

The observed asymmetry of the molecular disk can be 
partly attributed to the incomplete velocity coverage of the 
CO(2-l) emission in our data. As illustrated in Fig. 3, which 
shows the CO(2-l) spectrum observed towards the AGN, any 
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Fig. 4. a) V-H color map (gray scale in STmag magnitude system) of 
3C 236, with the integrated 12 CO(2 - 1) emission overlaid. Contour 
levels as in Fig. 2. b) Zoom into the central 3" of 3C 236. Contours 
represent the 12 CO(2 - 1) line maps of 3C 236 integrated in two chan- 
nels: blue (t£°-330 km s _1 <y<i)™) and red (i^<iKU™+485 km s -1 ). 
The contour levels are 3<r to 4.5cr in steps of 0.5 <x (cr=0. 13 Jy beam -1 
km s" 1 ) for the blue channel, and 3<x to la in steps of 0.5 a (<r=0. 13 Jy 
beam -1 km s -1 ) for the red channel. The dashed line shows the major 
axis of the molecular gas disk. (Aa, A<5)-offsets in arcsec are relative to 
the location of the AGN. The gray ellipses show the 1 mm beam. Color 
version available in electronic format. 



potential emission at velocities u-u„<-700 km s lies out- 
side the range covered by the PdBI narrow band correlator, 
which was initially centered in our observations around z = 
0.10054 (y (HEL)=30 129 km s 1 ). With this caveat in mind, 
the Gaussian fit to the CO spectrum of Fig. 3 provides an up- 
per limit to the systemic velocity (u sys ) of the source, assum- 
ing the plausible scenario where gas motions are driven by ro- 
tation at these radii. The shape of the CO line profile, which 
shows declining intensities at the low ('blue') velocity end of 
the spectrum, nevertheless indicates that the true i> sys is close 
to the adopted value: v-v a <-335 km s _1 , which corresponds to 
t£°(HEL)=29761+40 km s" 1 . The implied upper limit on the 
redshift z C o=0.09927+0.0002 is significantly smaller than the 
redshift adopted by Hill et al. (1996), who based their estimate 



on Paa and Ha measurements (z = 0.10054, Hill et al. 1996, 
and NED) 3 . The SDSS has a more recent optical spectrum of 
3C 236 available. The redshift of the source based on the SDSS 
spectrum is z = 0.0991+0.0001 (y sys (HEL)=29 740+20 km s -1 ), 
which is also consistent within the errors with zco- In the light of 
the new redshift value derived from CO, we analyze in Section 5 
the evidence of outflow signatures in different tracers of the ISM 
of3C236. 

Figure 2 shows a CO emission component SE of the AGN, 
unrelated to the rest of the CO. The size of this component is 
larger than the beam, and the flux is above the 3<x levels, an 
indication that this feature could be real. This component falls 
outside the outer dust disk and beyond the radio emission of the 
CSS source. It is also misaligned with the radio jets. There are 
no features in the VLBI, optical and UV maps at the location of 
this component, so its origin is unknown. 

Figure 4a shows the 12 CO(2 - 1) map superposed on the V- 
H color map made from HST imaging. The color map shows two 
prominent dust lane features that extend over significantly differ- 
ent spatial scales: the outer dust lane visualizes a highly-inclined 
gas disk of ~4.5" (8.1 kpc)-major axis diameter oriented along 
PA^55°. The inner dust lane pictures a much smaller disk-like 
structure of ~1" (1.8 kpc)-major axis diameter oriented along 
PA^30°. The two concentric disks are thus misaligned, a char- 
acteristic already noted by O'Dea et al. (2001) and interpreted as 
a signature of discrete accretion events. The inner disk is roughly 
perpendicular to the inner radio jet, oriented along PA^302°. 

4.2. Kinematics of molecular gas 

As shown in Fig. 4a, the CO disk seems to be closely linked 
to the inner dusty disk: to the limit of our sensitivity, CO 
emission is not detected in the outer disk. To better illustrate 
this association, we show in Fig. 4b the CO emission inte- 
grated in two channels defined to cover the ranges correspond- 
ing to blue velocities (u|r y °-330 km s -1 <y<i&°) and red veloci- 
ties (^ y °<u<if y °+485 km s -1 ). The velocity structure of the CO 
disk is spatially resolved in the maps: the emission peaks in the 
red and blue lobes are separated by ^0.45", i.e., a significant 
fraction of the PdBI beam. Lower-level yet significant emission, 
which stands out more clearly in Fig. 4b than in Fig. 4a thanks 
to the narrower velocity channels used to generate the maps, ex- 
tends up to radial distances ~0.7" (1.3 kpc) from the AGN. The 
line joining the peaks of the blue and red lobes has a PA=28.07°, 
which is consistent with the orientation of the inner dusty disk. 

Figure 5 shows the position-velocity (p-v) diagram taken 
along this line (PA=28.07°), identified as the kinematical major 
axis of the CO disk. The kinematical pattern of the p-v diagram 
reveals the signature of a spatially-resolved rotating molecular 
gas disk. At the NW edge of the CO disk, the radial velocity 
reaches ~400 km s -1 . If we assume that the inclination of the gas 
disk is 60°, based on the orientation of the radio source (Schilizzi 
et al. 2001), we derive that u, ot ~460 km s -1 at R~\.2 kpc. This is 
similar to the typical range of rotation velocities derived in other 
early-type galaxies (FWZI<800 km s ; e.g., Rrips et al. 2010; 
Crocker et al. 2012, and references therein). 

Based on a spherical mass distribution model, the dynam- 
ical mass Md yn inside R, derived as Md yn -Rxv^ ot /G (where G 
is the gravitational constant, R is the de-projected radius of the 



3 Hill et al. (1996) already noted that their redshift determination was 
significantly different from the previous estimate of Sandage (1967), 
who measured z = 0.0989 ± 0.0001, i.e., in close agreement with the 
CO-based value. 
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Fig. 5. Position-velocity diagram of the 12 CO(2 - 1) emission along the 
major axis of the molecular gas disk of 3C 236 (dashed line in Fig. 
4). Positions (Ax) are relative to the AGN. Velocities are relative to the 
CO u™. The contour levels are -2<x, 2<x to 4.5<x in steps of 0.5<x, with 
<r=3 mjy beam 1 km s -1 . 



disk, and u rot is the de-projected radial velocity at the edge of 
the disk), is ~6.1xl0 10 Mq. Accounting for the mass of molec- 
ular gas derived in Sect 4.3, we determine that the spheroidal 
stellar mass (M sph ) is ~5.8xl0 10 Mq inside R-l.2 kpc if the 
contribution from dark matter to Md yn is neglected. Not surpris- 
ingly, this value falls short of accounting for the total stellar mass 
M» of 3C 236, which amounts to 9.2 x 10 11 Mq according to 
the stellar population models of Tadhunter et al. (2011). This 
value, together with the available estimates of the supermassive 
black hole mass (Mbh) can be used to evaluate if 3C 236 fits 
within the known Mbh-M sp h relation followed by different galaxy 
populations. The two estimates for Mbh) in 3C 236 published 
by Marchesini et al. (2004) (M b h=4.2xl0 8 Mq) and Mezcua 
et al. (201 1) (M b h=3.2xl0 8 Mq) are in rough agreement. Taking 
the average of both estimates, the implied Mbh/M sp h ratio in 
3C 236 is -0.05%. This is compatible within the errors with the 
Mbh/M sp h ratio predicted for radio-loud AGNs at z~0.1, based 
on to the redshift-dependent law found by McLure et al. (2006) 
(see Eq. 3 of their paper). 



Strong emission from mid-IR H2 rotational lines, a tracer of 
warm (7\>100 K) and diffuse (n(H 2 )~10 3 cm- 3 ) molecular hy- 
drogen, has also been detected in 3C 236. The galaxy qualifies 
as a MOlecular Hydrogen Emission Galaxy (MOHEG) due to its 
large H 2 to PAH(7.7/mi) luminosity ratio L(H 2 )/L(7.7yum)>0.33, 
i.e., well beyond the limit defined by Ogle et al. (2010): 
L(H 2 )/L(7.7 y um)>0.04. Dasyra & Combes (2011) used Spitzer 
IRS to study the warm-H 2 emission lines ((0-0) SO 28.22 //m, 
(0-0) SI 17.04 fim, (0-0) S2 12.28 /mi, (0-0) S3 9.66 //m) of 
3C 236. Based on the flux of the SI, S2 and S3 lines, Dasyra 
& Combes (2011) computed that the total mass of warm-H 2 in 
3C 236 is 6.10xl0 7 Mq, with a single excitation temperature 
of T=345 K. Using the same data set, Guillard et al. (2012) 
fitted a LTE model of three excitation temperature components 
to the emission of the warm-H 2 in 3C 236. According to their 
model, the warm-H 2 phase is composed of (1.6 ± 0.5)xl0 9 Mq 
of gas at 7=100 K, plus (9.5 + 4.0)xl0 7 Mq at 7=241 K and 
7.3 ± 1.8 x 10 5 Mq for 7=1045 K. The warm-H 2 phase is more 
likely to have a distribution of temperatures, rather than one fixed 
temperature. Hence, we adopt the Guillard et al. (2012) model in 
our estimate of the warm-to-cold H 2 ratio in 3C 236, for which 
we obtain a value of 0.81. Gas-rich star forming galaxies usually 
present a warm-to-cold H 2 ratio of 0.01-0.1 (e.g., Higdon et al. 
2006; Roussel et al. 2007). Ogle et al. (2010) show that this ra- 
tio is significantly higher in powerful radio galaxies classified 
as MOHEGs, where the ratio ranges from -0.2 to ~2. The ratio 
measured in 3C 236 lies within the values found in MOHEGs. 

Although the optical line ratios measured in 3C 236 are con- 
sistent with those of low excitation radio AGNs (LERAGN, 
Smolcic 2009; Hardcastle et al. 2006), the M(H 2 )/M, ratio 
(-0.23%) and the M b h/M sph ratio (-0.05%) measured in 3C 236 
are similar to the values found in high excitation radio AGNs 
(HERAGN, Smolcic & Riechers 2011). The mixed properties 
of 3C 236 could be a consequence of the recent merger-driven 
reactivation of star formation and AGN activities in this source 4 . 



5. Evidence of outflow in 3C 236 

Based on the HI spectrum of 3C 236 obtained with the WSRT, 
which shows a 2000 km s _1 -wide absorption profile, Morganti 
et al. (2005) identified in this source an extreme outflow of 
atomic hydrogen. In the light of the new redshift of the source 
determined from CO, we discuss below the existence of outflow 
signatures that can be found from the analysis of the kinematics 
of different ISM components of 3C 236. 



4.3. Mass of molecular gas 

The total flux obtained within the 21" primary beam field of 
the PdBI is 4 Jy km s _1 . Assuming a brightness temperature ra- 
tio CO(2-1)/CO(1-0)-1, we can use the equations of Solomon 
et al. (1997) to derive the equivalent CO( 1-0) luminosity: L'co = 
4.6xl0 8 Kkm s _1 pc~ 2 . If we apply the Galactic ratio of H 2 -mass 
to CO-luminosity (4.6 Mq/K km s _1 pc 2 , Solomon et al. 1987), 
we estimate that the H 2 mass in 3C 236 is M(H 2 )=2.1xl0 9 Mq. 
Including the mass of helium, the corresponding total molecu- 
lar gas mass is 2.9 x 10 9 Mq. The average column density of 
hydrogen atoms derived from CO is /Vh.co = 9.4xl0 22 cm -2 . 
The Chandra spectrum of 3C 236, which has a spatial resolution 
(-0.6") similar to our observations, is best fitted by a power- 
law plus an absorption component with /Vh.co = 2xl0 22 cm -2 
(Birkinshaw et al., 2012 in prep.). The column densities derived 
from CO and X-rays are thus consistent within a factor 4-5. 



5. 1 . Cold molecular gas: CO line emission 

The kinematics of the spatially-resolved CO disk, analyzed in 
Sect. 4, can be explained by circular rotation around the AGN. 
Our data do not sample the extreme blueshifted end of the H i 
outflow (beyond v-v^®<-332 km s~'). We nevertheless note that 
the data show no signs of significant emission at the high end 
represented by extreme red velocities (883 km s~ 1 >i>-i;!r y °>368 
km s~'), which are fully covered in our observations. The im- 
plied upper limit to the molecular mass of the red lobe of the 
outflow is <4.4xl0 8 Mq. We used the properties of the molecu- 
lar outflow of Mrk 23 1 as an extreme upper limit to what could 
be expected in 3C 236: if we adopt in 3C 236 the same scal- 
ing for the outflow/disk mass ratio measured in Mrk 231 (i.e., 

4 LERAGN are usually hosted by galaxies on their final stages of the 
mass assembly, with older stellar populations and redder colors com- 
pared to their HERAGN counterparts. 
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Fig. 6. Hi absorption (Morganti et al. 2005), H2 S(l) 17 fim emission 
(Dasyra & Combes 201 1), and Gaussian fit of the 12 CO(2 - 1) emission 
of 3C 236. For comparison, the spectra have been normalized, and the 
Hi absorption is shown inverted. The red vertical line marks the u™. 
Top panel: Comparison of the CO emission (red, dashed line) and Hi 
absorption. The black line shows the integrated H 1 absorption spectrum 
of 3C 236. The three components of the fit to the H 1 spectrum are rep- 
resented with blue lines. Bottom panel: Comparison of the CO emis- 
sion (red, dashed line), H 2 S(l) 17 fim emission (green line), and the 
H 1 intermediate width component (blue line). The horizontal green line 
shows the FWHM of the H2 S(l) 17 jim line, corrected for instrumental 
broadening. 



M(H 2 ) outflow /M(H 2 ) disk ~0.06, Feruglio et al. 2010; Cicone et al. 
2012), the expected H 2 mass in the outflow would amount to 
1.2x10 s Mq. This is a still factor ^4 below the detection limit 
derived above for 3C 236. We can therefore conclude that the 
presence of a molecular outflow that would share the extreme 
properties measured in Mrk 23 1 would have remained unnoticed 
in 3C 236 by our observations. 

5.2. Atomic gas: HI line absorption 

Figure 6 shows the integrated H 1 absorption profile of 3C 236 
(Morganti et al. 2005). The spectrum has been inverted for the 
sake of comparison with the CO emission profile. The H 1 line 
profile shows three distinct velocity components hereafter re- 
ferred to as narrow, intermediate, and outflow components. The 
narrow and intermediate components were already detected by 
van Gorkom et al. (1989) using the VLA. The outflow com- 
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Fig. 7. Optical rest-frame spectrum of 3C 236 from SDSS 
(zsdss=0.0991). For clarity, we only show the region of the Hy8 
and [O m]/l/14959,5007 A emission lines. 

ponent has been confirmed with EVLA observations (A array, 
Morganti et al. in prep.). The velocity centroids and widths de- 
rived from the Gaussian fits to the narrow and intermediate com- 
ponents {narrow: u(HEL)=29 828 km s" 1 , FWHM-80 km s" 1 ; 
intermediate: u(HEL)=29 846 km s"\ FWHM-300 km s 1 ) dif- 
fer significantly from the parameters of the outflow component 
(u(HEL)=29474km s ~\ FWHM-l 100 km s 1 ). As argued be- 
low, the H i spectrum can be explained by two kinematically dis- 
tinct systems: a rotating disk and an outflow. 

The main limitation of Morganti et al. (2005) data, if we 
are to pinpoint the origin of the different components of the 
Hi spectrum, is their lack of spatial resolution. The compari- 
son between the integrated H i and CO profiles can be used to 
discuss the possible location of the H i absorbers, however. As 
shown in Fig. 6, the velocity range covered by the narrow and 
intermediate components is similar to the one covered by CO 
emission. This similarity can be taken as indirect evidence that 
the deep H i absorption is explained by a rotating H i structure 
that is likely concomitant with the molecular disk. This would 
leave out the evidence of H i outflow only to the most extreme 
velocities (-1000 km s 1 <y-yf°<-500 km s 1 ) 5 . The new VLBI 
map of 3C 236 presented by (Struve & Conway 2012) has re- 
cently confirmed this scenario by showing that the H i absorption 
of the narrow and intermediate component come from a rotat- 
ing structure characterized by an orientation almost identical to 
that of the CO disk. The location of the Hi outflow remains 
unknown, however, because the velocity range of the outflow 
lies well beyond the bandwidth covered by the data of Struve & 
Conway (2012). 

5.3. Warm molecular gas: H 2 line emission 

Figure 6 (bottom panel) shows the H 2 S(l) emission line profile 
detected in 3C 236 by Dasyra & Combes (201 1) 6 . A comparison 
of the S(l) line with the Gaussian fit to the 12 CO(2 - 1) spectrum 
and with the H i intermediate component indicates that the line 
profiles of these tracers are similar (Fig. 6, bottom panel). The 

5 We note that previous claims of an extreme H i outflow in 3C 236 
were mainly due to the wrong value assumed for u sys . 

6 The S(l), S(2) and S(3) lines detected in 3C 236 show similar 
FWHM (Dasyra & Combes 2011). 
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Fig. 8. Zooms of the H/3 (top panel) and [O m]/15007 A (bottom panel) 
emission lines (black, solid line), with the 12 CO(2 - 1) emission (red, 
dashed line), and Hi outflow absorption (blue, dotted line) overlaid. 
For comparison, all fluxes have been normalized, and the H i absorp- 
tion has been inverted. The horizontal dashed line shows the FWHM 
of I2 CO(2 - 1) emission. The vertical dashed line marks the u™. The 
shaded areas mark the emission of the optical lines beyond ±400 km 
s" 1 , above the continuum (yellow) and 3cr (green) levels. Color version 
available in electronic format. 



velocity centroids and deconvolved widths of the lines are al- 
most identical within the errors (Au~50 and 60 km s _I for H i and 
H2 centroids respectively). Even though the apparent FWHM of 
the SI line (750+75 km s~') is ~50% larger compared to the 
CO line and the H 1 intermediate component, this difference dis- 
appears if the instrumental broadening of the S(l) line is taken 
into account (FWHM=582±60 km s _1 , corrected for instrumen- 
tal broadening). This similarity can be taken as indirect evidence 
that the warm-H2 emission of 3C 236 is probably generated in 
the same rotating disk as the CO emission. 



'blue' velocities: |u-yf°|>400 km s _1 , i.e., well beyond the range 
allowed by the rotation of the disk (see discussion in Sect. 4). 
While the velocity centroids of the optical lines are similar to 
ifo, their widths (FWHM-900 km s 1 , FWZI-2000 km s 1 ) 
are ~3 times the value measured for the 12 CO(2 - 1) line. 

The imprint of the outflow is clearly detected on the red 
and the blue wings of the optical lines, suggesting that, unlike 
the emission from the star-forming disk, the emission from the 
outflowing gas is not as heavily extincted (see discussion in 
Sect. 6.1.1). The similar FWZI of the line wings suggests that 
the red wing is produced in the receding side of the same outflow 
system. The lack of spatial resolution prevents us from clearly 
identifying the origin of the outflow, however. A comparison 
with the H 1 outflow component shows that the blue wing of the 
ionized gas emission lines covers a comparable velocity range, 
which suggests that the ionized gas outflow may have started to 
recombine and form H 1 (e.g., Morganti et al. 2003b). 



6. Star-formation properties of 3C 236 

Nesvadba et al. (2010) studied the relation between the SFR sur- 
face density, and the molecular gas surface density derived from 
CO in a sample of radio galaxies. The SFR was estimated based 
on the 7.7 fim PAH feature emission, except for 3C 326 N, where 
the 70 fim continuum emission was also used as a SFR tracer 
(Ogle et al. 2007). With all the caveats discussed by Nesvadba 
et al. (2010) in mind, an inspection of the Kennicut-Schmidt 
(KS) law shown in Fig. 9 of their paper indicates that powerful 
radio sources as systematically associated with PAH intensities 
that are a factor ~ 10-50 lower than normal galaxies for a given 
mass. This offset is suggestive of a lower star-formation effi- 
ciency {SFE=SFR/M(H 2 )) in MOHEGs compared to normal 
galaxies. 

Nesvadba et al. (2010) pointed to the influence of large-scale 
shocks as the agent responsible for increasing the turbulence of 
molecular gas (see also Nesvadba et al. 2011). This energy in- 
jection in the medium could inhibit star formation to a large ex- 
tent. In Sect 6.2, we study the SFE and the location of 3C 236 
in the KS diagram, and discuss the differences with respect to 
the results obtained by Nesvadba et al. (2010) in their sample 
of MOHEGs. To this aim we have used several tracers of the 
SFR available for 3C 236 (discussed in Sect 6.1). This is re- 
quired if we are to constrain the potential biases inherent to the 
different SFR calibrations. Some SFR tracers are sensitive to the 
presence of AGN radiation and jet-induced shocks in the ISM, 
as these may increase the flux of ionized gas emission lines, 
vary the shape of the continuum, and destroy ISM molecules 
like the PAH. A possible evolutionary scenario explaining the 
differences between the samples is discussed in Sect. 6.3. 



5.4. Ionized gas: optical line emission 

Figure 7 shows the optical emission lines H^ and [Ora] 
/f/14959,5007 A of 3C 236, taken from the SDSS archive. The 
three lines show a remarkable large width, also seen in the op- 
tical emission of Ha, [Nn]^6548,6584A, [S n]AA67 13,6731 A 
lines, as well as in the IR lines [Nen] and [Nem] (Dasyra & 
Combes 2011; Guillard et al. 2012). Figure 8 shows a zoom on 
the H/3 and [Oin] /15007 A lines, with the Gaussian fit to the 
12 CO(2 - 1) spectrum and the Hi outflow component overlaid. 
This figure shows that a sizable fraction of the total emission 
of the lines (above 3<x-levels) is detected at extreme 'red' and 



6.1. Star-formation rate estimations 
6.1 .1 . Optical/Ultraviolet photometry 

Tremblay et al. (2010) and O'Dea et al. (2001) used HST high 
resolution imaging in the FUV, NUV, U, and V bands to study 
the star formation in 3C 236. As illustrated in Fig. 9, the FUV 
emission consists of five knots distributed in the inner 4.5 kpc of 
3C 236: there is one central knot located at the AGN and four 
knots arranged in an arc along the outer edge of the inner dusty 
disk described in Sect. 4. Tremblay et al. (2010) concluded that 
all the knots are associated with star-formation episodes yet of 
significantly different ages: ~10 9 yr for the nucleus, and ~10 7 yr 



8 



A. Labiano et al.: Fueling the central engine of radio galaxies 



1 -1 1 




O 






C) ■ ■ 


"■ *, 




■' 








n ■ n 


. o - 


, / 1 







Aa(arcsec) 




L • - 4) ■- 'I 



1 0.5 -0.5 -1 

Aoc(arcsec) 

Fig. 9. Top panel: HST FUV image of the star forming regions of 
3C 236, with the integrated 12 CO(2 - 1) emission overlaid. Bottom 
panel: zoom into the central 3" of 3C 236. The contours correspond 
to the 12 CO(2 - 1) emission integrated in the red (thin, black line) and 
blue (thick, gray line) channels. Contour levels as in Fig. 4. Color ver- 
sion available in electronic format. HST image in counts s" 1 . 



for the arc. Using Starburst99 stellar population synthesis mod- 
els to fit the data, Tremblay et al. (2010) determined that the 
UV-derived SFR of 3C 236 is 7.51 Mq yr" 1 . 

Figure 9 shows that the star forming regions of 3C 236 lie 
on the nearest edge of the CO and inner dust disk, while the rest 
of the galaxy shows no UV emission (except for the compact 
nucleus). This morphology, strongly asymmetrical with respect 
to the galaxy major axis, suggests that there is a strong extinc- 
tion in the inner dust disk. If we calculate the extinction using 
the N(H) values derived from CO and X-ray data (Sect. 4.3), 
the Bohlin et al. (1978) relationships yield Ay=10 and A v =50 
respectively. Therefore, the UV images of 3C 236 are probably 
missing UV emission from the star forming regions embedded 
in the inner disk and, thus, the SFR estimate of Tremblay et al. 



(2010) is probably a lower limit of the total SFR of 3C 236. The 
emission detected in the UV knots observed would arise from 
the outermost, low-extinction, layers of the inner dust disk. 

The extinction of the inner dust disk does not seem to sig- 
nificantly absorb the emission of the ionized gas outflow, which 
shows both red and blue wings in the spectrum (Sect. 5.4). Based 
on the emission line fluxes from the SDSS spectrum of 3C 236, 
the Cardelli et al. (1989) extinction law yields significantly lower 
extinction values: Ay =1.4 (using E(B - V)=0.456, Tremblay 
et al. 2010). This suggests that the outflow is formed in a dif- 
ferent region of 3C 236, probably the compact nucleus, which is 
visible in the UV images. High-resolution spectroscopy would 
be nevertheless needed to identify the location of the outflow. 

UV and optical emission of powerful radio galaxies may 
have a significant contribution from the AGN. This may cause 
a severe overestimation of the SFR if this contribution is not de- 
ducted. Tremblay et al. (2010) discussed that the contribution 
from emission lines to their UV and optical bands is negligi- 
ble for 3C 236, however. A similar conclusion was reached by 
Holt et al. (2007), who fitted the shape of the continuum emis- 
sion subjacent to the optical spectrum of 3C 236 inside the range 
3000-7000 A. In their analysis, Holt et al. (2007) found no evi- 
dence of a significant AGN component. 

Some CSS sources are known to show UV nebulosities that 
are roughly co-spatial with the region occupied by the expand- 
ing radio jet (Labiano et al. 2008). This association suggests that 
both emission features are related in these sources. In the case 
of 3C 236 however, the radio jets are perpendicular to the star- 
forming disk revealed by the HST images. This different geom- 
etry indicates that the UV and optical knots are mostly unrelated 
to the jet lobes in this source. We can therefore conclude that 
the UV and optical emission from the star-forming regions of 
3C 236 has a negligible contribution from the AGN. The value 
derived for the SFR (7.51 Mq yr -1 ) from these tracers is in all 
likelihood not significantly biased. 

6.1.2. Optical spectroscopy 

The Ha emission of a galaxy is indicative of its SFR (Kennicutt 
1998). Using the Ha line flux from SDSS (F Ha = 5.62 x 10~ 15 
erg s _1 cm 2 , Buttiglione et al. 2009), and the Ha-SFR relation- 
ships of Kennicutt (1998), we obtain a SFR=l.l Mq yr 1 for 
3C 236. If a galaxy harbors an AGN, the nuclear activity can 
alter the Ha flux through photoionization and/or shocks. To esti- 
mate the contributions from the AGN to the Ha emission, we use 
the so called diagnostic diagrams or Baldwin-Phillips-Terlevich 
(BPT) diagrams (Baldwin et al. 1981; Veilleux & Osterbrock 
1987). Buttiglione et al. (2009) measured the de-reddened ratios 
of the optical emission lines in the SDSS spectrum of 3C 236. 
A comparison of the values of [Oin]/Hy6, [Nn]/Ha, [Sn]/Ha 
and [Oi]/Ha from Buttiglione et al. (2009) with the results 
of Kewley et al. (2001) and Brinchmann et al. (2004) shows 
that the contribution of the AGN to the Ha flux of 3C 236 
is <10%. Removing 10% of the flux of the Ha line yields 
SFR— 1.0 Mq yr. Based on the same BPT diagrams, we find that 
the optical line ratios of 3C 236 are also compatible with ion- 
ization from shocks with velocities <300 km s . These shocks 
could be generated either by the AGN or by star formation pro- 
cesses. 

The 4000 A break (D4000, Balogh et al. 1999) is indicative 
of the SFR per unit of stellar mass (M») in AGN (Brinchmann 
et al. 2004). For 3C 236, D4000^1.8 (Holt et al. 2007), which 
yields S FR/M t ~lxlQ- n -lQ- 12 yr 1 . Using the total stellar 



9 



A. Labiano et al.: Fueling the central engine of radio galaxies 



mass of 3C 236 estimated by Tadhunter et al. (201 1), M«=10 12 
Mq, we obtain SFR-1-10 Mq yr -1 , consistent with the SFR 
estimations of Sect. 6.1.1. 

6.1.3. Mid-infrared continuum 

Calzetti et al. (2007) presented two SFR calibrations for their 
sample of nearby galaxies, one using 24 fim luminosity of the 
source, the other one combining the 24 //m and Ha luminosities 
(see also Kennicutt et al. 2009). Applying these calibrations to 
the 24 /an emission CF24^m=2.16xlO~ 12 erg s _1 cirT 2 , Dicken 
et al. 2010), and the Ha emission (Sect. 6.1.2) of 3C 236, we 
obtain S Ftf=6.62M yr" 1 and S FR=9 .63 ± 1 .7 Mq yr 1 respec- 
tively. 

Dicken et al. (2010) also measured the 70 /im flux of 3C 236 
with Spitzer. F 7o //m =2.77xlO~ 12 erg s _1 cm~ 2 . Based on the 70 
/im emission, we can establish a low and an upper limit (Dicken 
2011, private communication; see also Seymour et al. 2011) on 
the SFR of 3C 236: 3 Mq yv l <S FR<10 Mq yr 1 . 

The 24 /im and 70 //m emission can be increased by the heat- 
ing of dust by the AGN (Tadhunter et al. 2007), producing an 
overestimated SFR. However, Dicken et al. (2010) found that 
the 70 fim emission of 3C 236 is ~9 times larger than expected 
for a non-starburst galaxy with the same [O in] luminosity. They 
argue that this difference is attributable to star formation being 
the main contributor to the 70 //m flux. The SFR derived from 
the 24 pm and 70 //m are consistent with the SFR derived in 
Sect. 6.1.1, suggesting that the contribution from the AGN to 
these fluxes is small, and the SFR estimations from the IR con- 
tinuum are accurate. 

6.1.4. Mid-infrared spectroscopy 

Willett et al. (2010) used the [Nera]il5.6 //m and [Nen]^12.8 
fim fluxes to measure the SFR in a sample of Compact 
Symmetric Objects (see also Ho & Keto 2007). The Neon emis- 
sion lines of 3C 236 have fluxes F [NeII] = (0.92+0.06)xl0 10 
erg s" 1 cirT 2 and F [Nera] = (0.44+0.03)xl0 10 erg s _1 crrT 2 
(Spitzer spectrum, Guillard et al. 2012), which yield SFR = 
16±lM©yr _1 . 

To assess the AGN contribution to the total MIR luminosity 
of a galaxy, Willett et al. (2010) compare the [Oiv]/l25.9 yum / 
[Nen] /112.8 yum ratio with the equivalent width of the 6.2 yum 
PAH emission. We measured [Oiv]/[Nen]<0.2 and, EW 6 . 2 <170 
yum (3cr limit) from the Spitzer spectrum of 3C 236. Hence, the 
AGN contribution in 3C 236 is <10%. Genzel et al. (1998) used 
the 7.7 /im PAH strength, instead of the 6.2 yum, to study the 
contribution of the AGN to the galaxies in their sample. If we 
consider the 7.7 /im PAH strength of 3C 236, the contribution 
from the AGN is <20%, consistent with the samples of star- 
formation dominated galaxies (Genzel et al. 1998). Assuming 
that the AGN contributes in the same proportion to the MIR lu- 
minosity and the Neon lines, removing 20% of the Neon lines 
flux yields SFR=\2.9±Q.5 Mq yr -1 . 

Another indication of the AGN contribution is given by the 
ratio [Ne m]/[Ne n], which increases with the hardness of the ion- 
izing environment. For 3C 236, log [Neiii]/[Nen]=-0.3, similar 
to the mean value of the ratio in ULIRG and starburst galaxies 
(-0.35), and below the mean value for AGN (-0.07) and CSO 
(-0.16, Willett et al. 2010), suggesting a low contribution from 
the AGN in 3C 236. 

Pereira-Santaella et al. (2010) published a list of IR emis- 
sion line ratios for their sample of 426 active and HII galaxies. 



Based on their results, the [O iv]/[Ne n] and [Ne m]/[Ne n] ratios 
of 3C 236 are closer to the values observed in LINER-like and 
HII galaxies rather than in Seyferts or QSO. Therefore, the SFR 
of 3C 236, obtained with IR line ratios, is not largely overes- 
timated due to AGN effects. It is also consistent with the SFR 
values derived in Sect. 6.1.1. 



6.1.5. PAH emission 

The Spitzer spectrum of 3C 236 shows that the emis- 
sion of the PAH features is Fn /Jm =3.7xl0 -14 erg s _1 crrT 2 , 
F7.7^m<6.8xl0" 14 erg s _1 crrT 2 , and F 6 . 3/Jm <3.8xl0" 14 erg s _1 
cirT 2 (3cr limits, Dicken et al. 201 1; Guillard et al. 2012). 

Based on the 7.7 //m PAH emission limit, we used the 
Calzetti et al. (2007) equations (see also, Nesvadba et al. 2010) 
to establish a 3cr-limit on the SFR of 3C 236: SFR<0.23 Mq 
yr -1 . This upper limit falls below any other SFR estimation for 
3C 236, and is consistent with the SFR (which was also calcu- 
lated using the 7.7 /im emission) of the Nesvadba et al. (2010) 
radio galaxies with H2 masses ~10 9 Mq. 

Willett et al. (2010) used a second method to calculate the 
SFR in their sample, based on the correlation found between the 
Neon emission line luminosities and the 6.2 //m plus 11.3 //m 
PAH luminosities (see also Farrah et al. 2007). For 3C 236, the 
total flux of the 6.2 yum and 1 1.3 yum PAH emission is 3.7xl0~ 14 
erg s" 1 cm" 2 <F6.2 jU m+Fn.3 jU m<7.5xlO -14 erg s" 1 cm" 2 , yield- 
ing 1 1 Mq yr~' <SFR<22 Mq yr -1 . Leaving aside the estimate 
derived from the 7.7 yum PAH feature, we conclude that the av- 
erage value of the SFR obtained from all the different tracers 
discussed above is SFR~9.2 Mq yr _1 (the corresponding value 
from the 7.7 /im emission is thirty-five times lower). It is thus 
tempting to speculate that 7.7 ^m emission underestimates the 
SFR in 3C 236. 

O'Dowd et al. (2009) found that sources with an AGN com- 
ponent have weaker 7.7 /im emission than quiescent galaxies, 
consistent with the destruction of smaller PAH by shocks and/or 
radiation from the AGN. Modelization of the PAH interactions 
with the ISM show that PAH are destroyed by shocks with ve- 
locities > 125 km s _1 (Micelotta et al. 2010a,b, and references 
therein). The shock velocities measured in the warm-H2 phase of 
3C 236 are however too small (<30 km s"\ Guillard et al. 2012) 
to destroy the PAH 7 . Another mechanism of PAH destruction is 
by the AGN radiation field. Hard-UV and X-rays will destroy 
PAH at kiloparsec distances from the AGN (e.g., Siebenmorgen 
et al. 2004; Voit 1992). As discussed above, the star forming 
knots in 3C 236 are not found to be contaminated by the ef- 
fects of AGN radiation, however. Similarly, the PAH destruction 
from hard-UV and X-rays should in all likelihood affect only 
the nuclear region, leaving aside the star forming knots far from 
the AGN. The reasons behind the discrepancy between the SFR 
value derived from the 7.7 yum PAH and the rest of the tracers in 
3C 236 remain therefore unknown. 

6.2. The star-formation efficiency in 3C 236: 'normal' or 
'inhibited'? 

Figure 10 shows the SFR surface density (Zsfr =S FR/area) and 
C0M-H2 mass surface density (Z M m=M(H2)/area) for 3C 236 



7 PAH emission is likely produced in the PDR-like phase of molecu- 
lar gas at temperatures >100 K, which is also the major contributor to 
the bulk of H 2 line emission. 
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Table 2. SFR estimations of 3C 236. 



Region 


Tracer 


SFR 


M(H2) 


Area 


Refs. 






(Mq yr" 1 ) 


(10 s M Q ) 


(kpc z ) 




All 


Optical+ U V 


7.51 


21.3 


12.7 


1 


Unresolved 


Ho- 


1.0 


21.3 


12.7 


2,3 


Unresolved 


D400U 


1-10 


21.3 


12.7 


4,5 


Unresolved 


24 /im 


6.62 


21.3 


12.7 


6,7 


Unresolved 


24/jm+Ho' 


9.63+1.7 


21.3 


12.7 


2,6,7 


Unresolved 


70 /an 


3-10 


21.3 


12.7 


6,8 


Unresolved 


7.7 [im 


<U.z4 


O 1 'J 

li.j 


Iz. / 


"7 O 


Unresolved 


6.2+11.3 fim 


11-22 


21.3 


12.7 


10,11 


Unresolved 


;Neii]+[Nera] 


12.9+0.5 


21.3 


12.7 


10,11 


Region 


FUV 


6.19 


21.3 


12.7 


A 


Region 1 


FUV 


0.84 


<9.93 


3.0 


A 


Region 2 


FUV 


1.43 


5.49 


2.4 


A 


Region 3 


FUV 


1.27 


4.93 


2.5 


A 


Region 4 


FUV 


1.48 


3.66 


2.5 


A 


Region 5 


FUV 


1.17 


<12.5 


2.3 


A 


Limits are 3cr 


. Uncertainties 


in M(H 2 ) 


are -5%. 


SFR uncertain 



ties listed when available. 

References: 1-Tremblay et al. (2010), 2-Buttiglione et al. (2009), 
3-Kennicutt (1998), 4-Holt et al. (2007), 5-Brinchmann et al. (2004), 
6-Dicken et al. (2010), 7-Calzetti et al. (2007), 8-Seymour et al. (201 1), 
9-Guillard et al. (2012), 10-Dicken et al. (2011), 11-Willett et al. 
(2010), A- Appendix A of this work. 



and the Nesvadba et al. (2010) sample 8 . It also shows the canoni- 
cal KS-law fitted for normal star-forming galaxies (Roussel et al. 
2007; Kennicutt 1998). For the surface density calculations, we 
have assumed that all the star formation of 3C 236 is produced 
in a disk of radius 72=1.1" (deprojected area 12.7 kpc 2 ), which 
includes the molecular gas disk, and all the UV emission from 
the star-forming knots. Table 2 lists all the SFR estimations for 
3C 236, the areas considered, and the corresponding C0W-H2 
mass. 

Figure 10 shows that the Z S fr of 3C 236 is consistent with 
the Zsfr of normal star-forming galaxies at similar Smh2- In 
terms of SFE, 3C 236 shows SFEj,ci36-^ Gyr -1 . For a nor- 
mal galaxy with the same E M h2 as 3C 236, the KS-law yields 
£sfr=0.33 Mq yr -1 kpc -2 . Therefore, the implied S FE noima \=2 
Gyr" 1 . For the Nesvadba et al. (2010) sample, the SFE (0.02- 
0.04 Gyr -1 ) is 10-50 times lower than measured in normal 
galaxies. Compared to the powerful radio galaxies of Nesvadba 
et al. (2010), 3C 236 is an efficient star-forming radio galaxy. At 
the current rate, 3C 236 will deplete its H2 gas in fdepiete~2 Gyr. 
Our SFR calculations do not consider the extincted star-forming 
regions embedded in the inner dust disk (Sect. 6.1.1). The in- 
clusion of these regions in the SFR estimations would increase 
the SFE value of 3C 236 and thus the difference with the radio 
galaxies of Nesvadba et al. (2010). 

The HST images and CO map of 3C 236 provide enough res- 
olution to study the spatial distribution of the SFE. In Appendix 
A, we discuss how the Esfr and £ M H2 change along the molec- 
ular gas disk of 3C 236. Our analysis shows that all the star- 
forming regions of 3C 236 have a SFE consistent with the KS- 
law of normal galaxies. 3C 236 shows a high global SFE com- 
pared with other radio galaxies. This property cannot be at- 
tributed to any particular location in the disk. 



8 The data available from Nesvadba et al. (2010) do not in- 
clude warm-FL. Therefore, we considered only the cold-FL mass for 
the comparison of star-formation laws, which yields log E MH 2-8.2 
Mq kpc -2 . Adding the cold and warm-H 2 masses of 3C 236 yields log 
2 M H2^8.5 Mq kpc- 2 . 
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Fig. 10. Xsfr and S MH 2 of 3C 236 (triangles), the Nesvadba et al. 
(2010) sample (diamonds), and young sources 4C 12.50 and 4C 31.04 
(squares). The bottom panel is a zoom on the area marked with a dot- 
ted box, where most of the SFR estimations of 3C 236 fall. Solid line: 
best-fit of the KS-law from Kennicutt (1998). Dashed lines: dispersion 
around the KS-law best fit for normal star-forming galaxies (Roussel 
et al. 2007; Kennicutt 1998). 

6.3. AGN feedback in radio-galaxies: an evolutionary path? 

3C 236 has undergone two epochs of AGN activity, with a 
~ 1 7 yr phase of inactivity in-between. After this dormant phase, 
an accretion event (e.g., a minor merger 9 or re-settlement of 
the gas) increased the amount of gas, reactivated the star for- 
mation and triggered the (~10 5 yr old) CSS source (Tremblay 
et al. 2010; O'Dea et al. 2001). With the AGN inactive at the 
time of the accretion event, the lack of large-scale radio jets, and 
the lobes disconnected to the core, the feedback effects of the 
old radio source could not affect the newly acquired gas. If an 
outflow was created at the time of the CSS ignition, it would 
have needed to travel at an average velocity of 20 000 km s -1 to 
reach the star-forming knots seeing in the UV images (located 
at ~2 kpc from the nucleus) in 10 5 years. This speed is ~ 10-20 
times higher than the fastest outflows observed in radio galaxies. 
Then, such scenario can be excluded. The molecular gas in the 
knots has thus not been affected by the CSS source feedback. 
Therefore, the knots form stars at a normal rate. Within 10 6 — 10 7 
yr, the feedback effects may nevertheless reach the knots and in- 

9 O'Dea et al. (2001) detect a small companion galaxy at 10" from 
the nucleus of 3C 236. 
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hibit the star formation, lowering the SFE to the values measured 
for the Nesvadba et al. (2010) radio sources (which have ages of 
10 7 -10 8 years). 

In the scenario described above, where the effects of AGN 
feedback are closely related to the age of the radio source, 
we may expect to measure SFE values typical of normal star- 
forming galaxies in the hosts of other young radio sources. 
Willett et al. (2010) measured the SFR of 9 young sources 
using the 6.2 and 11.3 //m PAH emission (as shown in Sect. 
6.1.5 for 3C 236) and obtained values from 0.8 to 47.8 M© 
yr -1 . Two of these have published cold H2 mass estimations: 
4C 12.50 (SFR=31.5 M© yr" 1 , M(H 2 )=(1.0+0.3)xl0 10 M©, 
Dasyra & Combes 2012) and 4C 31.04 (SFR=6A M© yr 1 , 
M(H 2 )=0.5xl0 10 M Q , Garcia-Burillo et al. 2007). Figure 10 
shows that the SFE of both sources are comparable to the SFE of 
3C 236 and thus consistent with the KS-law of normal galaxies. 
Although the data are scarce, young radio galaxies seem to have 
higher SFE than old radio galaxies. 

Even though the detection of H 1 and ionized gas outflows 
confirms the presence of AGN feedback in the host galaxies of 
young radio sources, like the ones examined in this work, the 
kinematics of the cold molecular gas and the SFE of the hosts 
seem, however, mildly affected. It is possible that, due to the 
extreme youth of the AGN, the effects of feedback in young 
sources have not had the time to propagate to such large scales 
as in the evolved radio sources. 



7. Summary and conclusions 

We have used the IRAM PdBI to study the distribution and kine- 
matics of molecular gas of the nearby (z~0.1) FR II radio galaxy 
3C 236, by imaging with high spatial resolution (0.6") the emis- 
sion of the 2-1 line of 12 CO in the circumnuclear region of 
the galaxy. AGN activity has been recently (~10 5 yr ago) re- 
activated in 3C 236 , triggered by a merger episode (that oc- 
curred ~10 7 yr ago). Previous observations identified in 3C 236 
one of the most extreme H 1 outflows thus far discovered in a ra- 
dio galaxy. The new observations presented in this paper have 
been analyzed to search for the footprints of AGN feedback on 
the molecular ISM of 3C 236. In particular, we have looked 
for evidence of outflow motions in the molecular gas kinemat- 
ics. Furthermore, we have derived the SFE of molecular gas in 
3C 236 and compared this value with the SFE obtained in differ- 
ent populations of radio galaxies, including young and evolved 
radio sources. We have investigated if the star-formation proper- 
ties of 3C 236 deviate from the KS relation followed by normal 
star-forming galaxies. 

We summarize below the main results and conclusions of 
this work: 

- The CO emission comes from a spatially resolved 2.6 kpc- 
diameter disk that has a gas mass M(H2)=2.1xl0 9 M©. 
The CO disk is linked to the inner region of a highly- 
inclined dusty disk identified in the V-H HST color image 
of the galaxy. The molecular disk is fueling a star formation 
episode singled out by a chain of strong emission knots de- 
tected at UV wavelengths in the HST pictures of 3C 236 . 

- The kinematics of the CO disk are dominated by circu- 
lar rotation and show no indications of high-velocities at- 
tributable to outflows. However, based on the limits im- 
posed by the sensitivity and velocity coverage of our data, 
we cannot exclude in 3C 236 the existence of a molecu- 
lar gas outflow comparable to the one detected in Mrk 23 1 . 



The gas disk rotates regularly around the AGN core. The lat- 
ter is identified by prominent radio continuum emission de- 
tected at 3 mm and 1 mm wavelengths in our map. Based 
on the observed CO kinematics we determine an upper 
limit for the redshift of the source z C o=0.09927±0.0002 
(u!r°=29761+40 km s _1 ). The new value of u sys is signifi- 
cantly blue-shifted (>350 kms~') with respect to the value 
of the systemic velocity previously reported in the literature 
for the source. 

- In the light of the new redshift value, the bulk of the deep H 1 
absorption can be simply explained by a rotating H 1 struc- 
ture, leaving out the evidence of outflow only to the most 
extreme velocities (-1000 km s _1 <y-u^<-500 km s~'). As 
for the CO emitting gas, outflow signatures are also absent in 
the warm molecular gas emission traced by infrared H2 lines. 
A reanalysis of the ionized gas kinematics reveals neverthe- 
less the existence of a fast outflow visible in the high-velocity 
(~1000 km s~') blue and red-shifted line emission wings. 

- We have derived the SFR of 3C 236 using the different trac- 
ers available at optical, UV and MIR wavelengths. The most 
reliable SFR value comes from high-resolution UV-optical 
photometry, which in 3C 236 seems mostly unaffected by 
AGN contribution. Most of the other SFR tracers examined 
may have contributions from the AGN, however. The con- 
sistency of the SFR values obtained with all SFR tracers 
suggests nevertheless that the AGN contribution is kept low. 
High-resolution spectroscopy of the UV star-forming knots 
would be needed to accurately measure the AGN contribu- 
tion. 

- The SFE of the molecular disk, defined as S FE=S FR/M gas , 
is fully consistent with the value measured in normal galax- 
ies, i.e., systems that form stars following the canonical KS 
law. This result is at odds with the factor of 10-50 lower 
values that have been claimed to characterize the SFE of 
evolved powerful radio galaxies. As an explanation for the 
different SFE measured in young and evolved radio sources, 
the most likely scenario in 3C 236 suggests that none of the 
effects of AGN feedback have yet had the time to affect the 
kinematics or the SF properties of the molecular ISM due to 
the extreme youth of the AGN re-activation episode in this 
source. 

This simple evolutionary scenario remains to be confirmed 
using larger samples, and most notably, using different tracers of 
the SFR in radio galaxies. As a caveat, it should be pointed out 
that most of the estimates of the SFR in evolved radio sources 
have been thus far derived from the strength of the PAH bands 
at 7.7 /im. The low SFR values in these sources may be severely 
underestimated due to the destruction of PAH carriers by shocks 
and/or strong UV and X-ray fields; these ingredients are known 
to be in place particularly in radio-loud AGN. The reported dif- 
ferences in the S FE of young and evolved radio galaxies could 
simply reflect a mismatch in the calibration of the SFR tracers 
used. 
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Firstly, we smoothed the HST FUV image with the beam 
of the 12 CO(2 - 1) maps, to have the same resolution in both 
images. We then divided the disk in five regions, based on the 
UV emission and molecular gas content. Figure A. 1 shows the 
smoothed FUV image with the 12 CO(2 - 1) blue and red chan- 
nels overlaid, and the location of the five regions. Regions 1 
and 5 correspond to the most external UV emission with no 
I2 CO(2 - 1) detection above 3<x levels. Regions 2 and 4 have 
UV and 12 CO(2 - 1) emission. Region 3 corresponds to the cen- 
ter of the disk. We named "Region 0" the sum of regions 1 
through 5. 

SFR=7.51 Mq yr~' was calculated using all flux down to 
lcr over the background (aperture "All", Tremblay et al. 2010). 
To obtain the local SFR for each region, we divided 7.51 Mq 
yr~' by the flux in the "All" aperture, and multiplied by the flux 
of each region. The FU masses were measured for each region 
individually, using the 12 CO(2 - 1) emission map. The SFR and 
H2 contents for each region are listed in Table 2. The SFE of the 
star forming regions are S ,F£ri>0.8 Gyr -1 , SFEr2=3 Gyr -1 , 
SF£ R3 =3 Gyr 1 , SF£ R4 =4 Gyr -1 , and SFE R5 >0.9 Gyr" 1 , con- 
sistent with the SFE calculated using the "All" aperture, and 
the SFE of galaxies following the KS-law (e.g., SINGS sample, 
Kennicutt et al. 2003; Roussel et al. 2007). The SFE of 3C 236 
is not a local effect from one individual region, 3C 236 shows 
normal SFE all along its molecular gas disk. 



Appendix A: Spatially-resolved SFE of 3C 236 

The emission of the star forming regions and the CO of 3C 236 
are almost co-spatial, suggesting that 3C 236 is using the cold 
molecular gas traced by CO to form stars. The resolution of the 
HST images and the CO map are enough to study the SFE along 
the molecular gas disk of 3C 236. Therefore, we searched for the 
star-forming region responsible of the normal (KS-law) SFE of 
3C 236. 




Aa(arcsec) 



Fig. A.l. HST FUV image of the star forming regions of 3C 236, 
smoothed with the beam of CO, with the red and blue l2 CO(2 - 1) chan- 
nels overlaid. Contour levels as in Fig. 4. Rectangles show the schematic 
location of the different regions for the analysis of the SFR and H 2 dis- 
tribution. Color version available in electronic format. 
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